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Abstract 

Bose-Einstein correlations in relativistic heavy ion collisions are very different for opaque sources than 
for transparent ones. The Bose-Einstein radius parameters measured in two-particle correlation functions 
depend sensitively on the mean free path of the particles. In particular we find that the outward radius 
parameter for an opaque source is smaller than the sideward radius parameter for sufficiently short duration 
of emission. A long duration of emission can compensate the opacity reduction of the longitudinal radius 
parameter and explain the experimental measurements of very similar side- and outward radius parameters. 



I. INTRODUCTION 



Bose-Einstein interference of identical particles (pions, kaons, etc.), emitted from the collision zone in 
relativistic heavy ion collisions shows up in correlation functions and is an important tool for determining 
the source at freeze-out. It is commonly assumed that the source is cylindrical symmetric around the beam 
axis and transparent, meaning that the detector receives particles from all over the source. The radius 
parameter outwards or towards the detector, R Q , is found to be larger than the radius parameter, R s , 
sidewards or perpendicular to the detector 

Rl = R 2 S + (3 2 Jt 2 , (1) 

for particles of zero rapidity. The excess is due to the duration of emission, St, of the source in which 
particles with transverse momentum p± and outward velocity (3 — p±/m± travel a distance f3 St towards 
the detector. Distances perpendicular to this velocity like, for example, R s are not affected by the duration 
of emission and thus reflect the "true" transverse size of the source, i.e., the region over which the source 
can be considered homogeneous. When flow is included the relation (|l|) is still valid in the analyses of 
whereas strong flow coupled to surface emission can reduce R Q more than R s || as well as an ellipsoidal 
source 

Experimentally the HBT radius parameters, R s ,o,i,oi, are extracted from measurements of the two-pion 
and two-kaon correlation function by parametrizing it with the common gaussian form 

C(q s , q , m ) = 1 + Xcxp(-q 2 s Rl - q 2 R 2 - q 2 R 2 - 2q oqi R 2 ol ) . (2) 

Surprisingly, in relativistic heavy ion collisions the outward and sideward radius parameters are measured 
to be similar [P-|lCfl (and in a few cases the outward size is even measured to be smaller than the sideward 
size HQ] contradicting equation ([!])) within experimental uncertainty. According to equation (|l|) this implies 
that particles freeze-out suddenly, St <C Ri, as in a "flash" in particular when resonance life-times are 
included HQ. 

We want to point out that an opaque source emitting away from its surface naturally leads to R <C R s 
unless the duration of emission, St, is very long. There is therefore a strong correlation between the direction 
of the emitted particles and the emission zone - contrary to cylindrically symmetric transparent sources 
without flow. In our model the measured sideward radius parameter samples the full source size, while the 
measured outward radius parameter sample only a small surface region of the opaque source in the direction 
of the emitted particles. The opaque source is inspired by hydrodynamical as well as cascade calculations 
of particle emission in relativistic heavy ion collisions. In hydrodynamical calculations particles freeze-out 
at a hypersurface that generally does not move very much transversally until the very end of the freeze-out 
However, most hydrodynamical freeze-out mechanisms as well as the analysis in Q do not include the 
directional condition that particles can only be emitted away from the surface, though strong flow at the 
surface has a similar effect. In cascade codes the last interaction points are also found to be distributed in 
transverse direction around a mean value that does not change much with time fl^-|l6f| . The thickness of the 
emission layer is related to the particle mean free path X m f P , and is zero in hydrodynamical calculations but 
several fin's in cascade codes. We shall in section 2 consider the parameters R and A m / P as constants and 
in section 3 discuss effects of moving surfaces, temporally dependent emission layers and transverse flow. 

II. SURFACE EMISSION AND HBT RADIUS PARAMETERS 

To quantify our statements we calculate the radius parameters R s ,o,i for two simple sources commonly 



symmetric but longitudinally expanding source. However, we add the important requirement that the source 
should be opaque emitting away from a surface layer of thickness ~ A m / p . 

For the correlation function analysis of Bose-Einstcin interference from a source of size R we consider 
two particles emitted a distance apart with relative momentum q = (ki — k2) and average momentum, 
K = (ki + ka)/2. Typical heavy ion sources in nuclear collisions are of size R ~ 5 fm, so that interference 
occurs predominantly when q ~ h/R ~ 40 MeV/c. Since typical particle momenta are ki ~ K ~ 300 MeV, 
the particles escape almost parallel (see Figure), i.e., k\ ~ &2 — -K" 3* 9- The correlation function due to 
Bose-Einstein interference of identical particles from an incoherent source is (see, e.g., 0) 



Jd 4 x S(x, K) e"* 
/d 4 ! 5(i,K) 



c(q 5 K)=i + i j : rri: 1 2 , 0) 



where S"(x, K) is a function describing the phase space density of the emitting source. With qx = q ■ x — 
q-/3^' t, where /3fc=K/ is the pair velocity, one can, by expanding to second order in i?^ and compare 
to equation (|J), find the radius parameters Ri, i=s,o,l,ol, 

R\ = a{ Xl -fat), R 2 ol = ((x - 13 t) (xi - ft t)) , (4) 

where the source average of a quantity O and its fluctuation are defined by 

The reduction factor A in equation (||) may be due to long lived resonances coherence effects, incorrect 

Gamov corrections [[l7| or other effects. It is found to be A ~ 0.5 for pions and A ~ 0.9 for kaons. 
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FIG. 1. Cross section of the interaction region perpendicular to the longitudinal or z-direction. Particles have 
to penetrate a distance I out to the surface of the interaction region in order to escape and reach the detector. The 
bulk part of the emitted particles comes from a surface region of width Xmfp, the mean free path of the particle. 

Generally a particle is emitted and escapes after it has made its final interaction with any of the other 
particles in the source. If the source is dense it is more probable to emit from the surface than from the 
center. In Glauber theory particles are emitted from a surface layer with a depth of order of the mean free 
path, A m / P , with probability 

POO 

S(x, K) ~ exp(— / dx'an{x')) — exp(—l/X m f p ) . (6) 

J X 

Here the integral runs over the particle trajectory from emission point x to the detector. Assuming an 
average density n in the emission layer, the distance I, that the particle has to pass through is of order of 
the mean free path, X m f p = 1/an (see Fig. 1). A similar surface layer was introduced by Grassi et al. 
in order to improve the standard Cooper-Frye freeze-out in hydrodynamical calculations and to study the 
effect on transverse momentum spectra. 

In the following subsection we assume that the source size R and thickness of emission layer X m f p are 
constants in time. In section 3 will will generalize these results including temporal dependences in connection 
with the final freeze-out of the source. 
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A. Spherical Sources 



First we investigate a spherically symmetric source of radius R. Including the temporal evolution of the 
source emission function by St(t), the phase space density of the source will have the approximate form 

S(x, K) ~ 9(i? 2 - x 2 - y 2 - z 2 ) e- 1 '^ S t (t) . (7) 

For strict surface emissionF] (A m / P = 0) the source reduces to 

S(x,K) ~ cos 6» e(cos0)6(R 2 - x 2 - y 2 - z 2 ) S t (t) . (8) 

where 9 is the polar angle with respect to the outward direction along K; we will in the following orient our 
cartesian axes such that the x-axis is in the outward direction. The geometric factor cos 9 — x/i? suppresses 
the peripheral zones and the O(cos0) = 0(x) factor insures that particles are only emitted away from the 
surface, i.e., only particles from the surface layer of the half hemisphere directed towards the detector will 
reach it whereas particles from the other hemisphere will interact on their passage through the source. 

The temporal emission is determined by St(t). It is commonly approximated by a gaussian, St(t) ~ 
exp(— (t — io) 2 /2(5i 2 ), around the source mean life-time, tq with width or fluctuation, <5t, which is the 
duration of emission. These gaussian parameters approximate the average emission time, (t) and the variance 
or fluctuation, er(i), for a general source, respectively. 

Orienting our coordinate system with outward or x-axis along K we have @k=K-/Ek — (/3 o ,0,0) (see 
Figure 1) when the pair rapidity has been boosted longitudinally to their center-of-mass system, Y = 0. In 
that case I = {y^ R 2 — y 2 — z 2 — x). When X m fp <C R we obtain from equations (^j5|) 

1 29 

Rl = ^R 2 + ^A 2 m/p + P 2 a{t) , (9) 

R 2 S = ^R 2 - gA^/p' ( 10 ) 

}2 r>2 



Rf = Rt , (ii) 

and R i vanishes. Aside from the /3 2 a{t) term, we notice that an opaque (A m fp <C R) source has i? 2 = 2i? 2 /9, 
i.e., the outward radius parameter seems much smaller than the actual source size as well as R s . This is a 
simple geometrical effect that arises because no particles are observed from the dark side of the source (see 
Figure 1) leading to an expectation value of (x) that is displaced from the center and only slightly smaller 
than yj (x 2 ). For comparison, a transparent (A m / P = oo) and spherically symmetric source has the same 
extent in all directions and one finds R s = R a = Ri = R/yh. Including the term due to the duration of 
emission, /3 2 a(t), the outward HBT radius parameter is necessarily larger than the sideward. 

B. Cylindrical Sources 

At very high energies longitudinal expansion is important and we employ a locally thermal but longitudi- 
nally expanding (with Bjorken flow u = z/t) source as in 



* A spherical source similar to (]8|) with strict surface emission, Xmfp = 0, was applied to proton emission from 
excited nuclei [^9). As proton decay times are long only the finite duration of emission contribution was considered 
(i.e, the last term in equation (H).) 
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2, / m± cosh(Y* — rj) I 



S(x, K) ~ Q(I? -x'- y) exp ( ± ^ '± ) S T (r) . (12) 



Here, r = \/t 2 — z 2 is the invariant or proper time and r\ = 0.51n(t + z)/(t — z) the space-time rapidity. R 
is the transverse source size, m± the transverse mass and Y the rapidity of the particles. The distance that 
the particles have to pass through matter in the transverse direction is I ~ {\/ R 2 — y 2 — x )- Actually, when 
f] Y the particles pass a distance longitudinally, which will lead to a minor reduction of Ri, but as we 
concentrate on the transverse radius parameters we will ignore this effect in the following. From equations 
(HH) we °btain when \ m fp "C R 









)i? 2 + 


fi- 


-) 




32 ) 



K = \^-i^Y)R z + [^-^)K lfP + (%<r{T), (13) 



R l - -^R 2 - Q^mfp > ( 14 ) 
Rl * (r 2 ) — , (15) 

for Y = (which leads to i? ; = 0) and keeping only leading orders in <j(t) and T/m± (see, e.g., |lj for 
higher orders). Again we notice that for a short-lived (/3 G ct(t) 1 / 2 <C R) and opaque (X m fp ^ R) source 
the outward radius parameter seems much smaller than the actual source size as well as the sideward radius 
parameter, R Q — 0.22R — 0.39R S (see Figure 2). For a short-lived transparent source (A m / P = oo) the 
transverse radius parameters are the same R s — R = R/2. The experimental measurements find, that 
R ~ R s . This implies a very long duration of emission, f3 o-(t) 1 / 2 ~ R s , for an opaque source and is in 
sharp contrast to the very short duration of emission the same data would imply for a transparent source. 



III. TEMPORAL DEPENDENCE AND FREEZE-OUT 



Eventually all particles created in the high energy nuclear collisions must freeze-out. In cascade codes the 
emission changes from strict surface emission shortly after the collision starts to an emission layer, which 
increases in thickness to the size of the system. At this stage volume emission takes over and, aside from long 
lived resonance decays, the source freezes out |l^-|l^]. This picture is somewhat different in hydrodynamic 
codes, where the emission is strictly from the surface at all times, but the surface eventually moves inward and 
the source freeze out when the surface hits the center fl3| . However, just before freeze-out the hydrodynamic 
hypersurface moves inward with superluminous speed. Again this stage is better characterized by volume 
freeze-out rather than by surface emission. The hydrodynamic results may become more similar to those 
of the cascade models, if we include emission from a surface layer of finite thickness, as in fis|] , as well as a 
cascade among the hadronized particles. 
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FIG. 2. Sideward and outward HBT radius parameters as function of the relative thickness of surface layer 
thickness to source size. The small mean free path limits are given in equations (^) and ( p^| ) and for X m f p ^> R they 
approach R s — R = R/2. The duration of emission contribution 01 cr(r) is not included in R 2 ,. The strength of the 
opacity effect, AR 2 /R 2 , is estimated from equation (^) at the freeze-out time Tf. 



The freeze-out can be described by allowing the layer of emission or \ m fp to increase with time after the 
collision. According to the Bjorken scaling model the mean free path increases linearly with proper time r 

1 ttR 2 



•Wp(r) 



^mfp{Tf)- 



(16) 



ctti(t) adN/dy " " ' Tf 

since the particle density drops inversely with proper time due to the one-dimensional expansion in longi- 
tudinal direction. In cascade codes the size of the emission layer also increases with time up to freeze-out. 
Assuming such a linear increase a simple estimate of the HBT radius parameters can be made. From the 
definition in equation ^ we find, that the HBT radius parameters are the time-averaged values of the X m fp 
dependent HBT radius parameters 



R 2 = — 



Rj(\ mfp {r)) t dr = / R 2 (X mfp (r)) 



A m/p( T /) 



(17) 



'/ JO JO 

Thus the HBT radius parameters are the average values of the mean free path dependent radius parameters 
of Figure (2) up to X m f p (rf). If we include now the duration of emission the difference is 
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R a - R 




(R (X m fp) - R 2 s {\ m fp)) 



\2 



•mfp( T f) 



= -AR 2 +(3 2 a(r). 



(18) 



Here AR 2 is the average distance between R 2 S and R 2 in Figure 2, i.e., the area between the curves for R 2 
and R 2 up to A 2 n y p (i/) and divided by that same factor A 2 ^ (r/). The opacity difference AR 2 is plotted as 
dashed line in Figure (2). Though it decreases with increasing A m / P (T/), it is a substantial fraction of R 2 
for mean free paths at freeze-out of order the size of the system. 

It is instructive to consider the recent and relevant example of central Pb + Pb collisions at energy 160 
A-GeV. The NA44 data finds for the pion HBT radius parameters R s ~ R a ~ 4.5 - 5.0 fm and Ri ~ 5 — 6 
fm P]. The average transverse momentum was pj_ ~ 165 MeV such that (3 = p±/m± ~ 0.76. From the 
transverse momentum slopes of pions, kaon, protons and deuterium in |20| one finds a temperature T ~ 120 
MeV and transverse flow v ~ 0.6c. Larger transverse flow and smaller temperatures can, however, not be 
excluded. The transverse flow effect on the HBT radius parameters is small for these p± and v values |H] . 
If we assume that the source emits pions at a constant rate per surface element, we can estimate from the 
longitudinal HBT radius parameter, equation (|l5|), the freeze-out time Tf — (2 m^/T) 1 / 2 Ri <~ 11 fm/c. As 
seen from equations ( |l4| ) the sideward radius parameter is relatively less affected by the X^fv term and we 
extract an initial transverse source size R ~ 3 1 / 2 R s ~ 9 fm from the NA44 sideward HBT radius parameter. 

To determine the outward HBT radius parameter we have to estimate the mean free path, for example, 
from equation (fl6|). In central Pb + Pb collisions at 160 A-GeV the particle rapidity density is dN/dy ~ 600 
around midrapidity. If we take for the pion cross section av ~ 20 mb, the resulting mean free path is very 
small, A m / P (rj) ~ 1 fm. However, near the surface the density will be lower than the average density used in 
this estimate and we would rather expect A m / P (r/) ~ R since pion freeze-out occur when their mean free path 
is of order the size of the system. For example, AR 2 ~ 0.13 (0.07) R 2 for A m / P = 1.0 (2.0) R respectively and 
the opacity effect leads to a corresponding difference AR 2 ~ 16 (9) fm 2 . Equation ( |l8| ) clearly demonstrates 
the significance of the opacity effect in reducing the outward radius parameter with respect to the sideward. 
The duration of emission compensates this difference to some extent. Assuming constant particle emission 
per surface element up to freeze-out, we find (3 2 <j(t) = (3 2 t 2 /18 ~ 4 fm 2 . This contribution due to the 
duration of emission is thus smaller than the reduction due to the opacity contribution and consequently 
one would expect R < R s . However, a number of other fluctuations such as moving surfaces, short lived 
resonances, and other effects also add to R 2 

We notice that the transverse source size is only a few fm larger than the geometrical size of Pb ~ 7 fm. 
Some small expansion of the source seems to take place before final freeze-out. The average emission time is 
yj (t 2 ) — t//v2 ~ 8 fm/c. This result is entirely consistent with the freeze-out time needed to explain the 
enhancement in the ir~ /tt + ratio at low p± in the same Pb + Pb collisions due to Coulomb repulsion J22|. 



We have addressed the significant difference in HBT radius parameters between opaque and transparent 
sources and have calculated the radius parameters and their dependence on the mean free path for simple 
spherical and longitudinal expanding models. A short-lived and very opaque source or "black body" has 
much smaller outward than sideward radius parameter. However, a long duration of emission as well as 
other fluctuations can compensate the reduction due to the opacity effect leading to comparable R ~ R s 



IV. SUMMARY AND OUTLOOK 
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as measured in several experiments. In the case of central Pb + Pb collisions at 160 A-GeV the opacity 
reduction was estimated to be a significant fraction of the outward HBT radius parameter and larger than 
the contribution from a source emitting during all its life-time. Other fluctuations like moving surfaces, 
short lived resonances, and other effects do, however, also add to R% and we expect that this is the reason 
that the sideward and outward HBT radius parameters are measured to be very similar in relativistic heavy 
ion collisions. We find that pions do not appear as in a "flash" . In fact a long duration of emission of 
order the life-time of the source is possible and consistent with an outward HBT radius parameter smaller 
or comparable to the sideward due to the opacity effect. In contrast, a simple transparent source would 
necessarily have a very short duration of emission as implied by equation (1). 

The Cooper- Frye freeze-out condition in hydrodynamic models does not take the opacity effect into account 
and generally one finds considerably larger outward than sideward source radius parameters due to the long 
freeze-out time fl3|| . Cascade codes have implicitly opacities build in through rescatterings and do find a 
directional effect Jl^] , but also long duration of emission and mean free paths |ll| are found leading to larger 
outward radius parameters than sideward. 

Other effects like, e.g., transverse flow, can also reduce the outward HBT radius parameter more than the 
sideward. At larger transverse momenta the sideward and outward HBT radius parameters are reduced by 
transverse flow and the longitudinal scales by the factor l/m±. This may explain the decrease of all the 
HBT radius parameters with increasing p± as found in the NA44 experiments [|| . 

By studying the rapidity and p± dependence of the HBT radius parameters the duration of emission 
contribution, 1 cr(r), may be separated from the opacity effect, AR 2 , as well as other fluctuations. However, 
this may be non-trivial if the mean free path, duration of emission and resonance contributions are rapidity 
and p± dependent. Transverse flow may also add a strong p± dependence of the radius parameters but 
the magnitude of transverse flow can approximately be determined from transverse momentum spectra . 
HBT radius parameters and transverse flow increase for heavier ions colliding. The nuclear A-dependence 
can therefore provide additional information on freeze-out times, sizes, transverse flow and opacity of the 
sources. 
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